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The theory of impurities in excitonic insulator is investigated in the light of the recent experiments 
on hexaborides. First, we study the bound state around the impurity and find that the bound state 
emerges when ReA is positive. Second, we study the continuum state using Abrikosov-Gor'kov's 
^^ ' approach. We find that the energy gap is reduced strongly when ImA = 0. Finally, we solve 

rvj ' Bogoliubov-de Gennes equations for excitonic insulator numerically. We get the results consistent 

with the analytic ones. We also find that incomplete ferromagnetism appears in doped triplet 
excitonic insulator with impurity. We make a short qualitative discussion on the ferromagnetism of 
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r^ ^ doped hexaborides using our result. 
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(D ■ I. INTRODUCTION 



C/3 I The discovers of weak ferromagnetism in La-doped CaBg is one of the greatest surprises in the recent study of 

ferromagnetismEI. This material showsLferromagnetism with high Curie temperature about Tc ~ 600K and with small 
magnetic moment. Zhitopairsky et al.a claimed that the theory of the excitonic ferromagnetism, which was originally 
proposed by Volkov et aZHI, can account for this curious ferromagnetism. Their proposal has brought renewed interests 



cd 



r^ ' in excitonic insulator, which was vigorously studied theoretically in 1960s. j-, 

(— I After the proposal |by Zhitomirsky et at, some authors tried to improve their theory. Balents and Varmaa and 

Q Barzykin and Gor'ko'vu independently discussed that, to explain this ferromagnetism, the approximation of excitonic 
O ' insulator should be improved by considering the formation of superstructure or long-range Coulomb interaction. 

The effect of impurities will also be important for the excitonic insulator. In this paper, we study the effect of 
^-^ ^ impurities in excitonic insulators. [-■ 

^ ■ The destruction of excitonic insulator by the impurity was discussed by Zittarta3. He showed that the intraband 

^^ ' scattering by impurity reduces the energy gap just as the magnetic impurity in superconductor. In this paper we 
CT^ study the interband scattering by impurity. We show that the interband scattering cause the formation of the bound 
^T^ ' states and the reduction of the energy gap. We also show that this interband scattering breaks the U{1) symmetry 
^sij \ of hamiltonian. This is due to the fact that the phase of the order parameter is determined by the difference of phase 
f—v between conduction and valence bands. The reduction of the energy gap and the energy of the the bound state is 
f~*) controlled by the phase of order parameter. The phase dependence of the energy gap and the bound state brings the 
*>*; ferromagnetic ordering of doped electrons in triplet excitonic insulator, because the phase of order parameter differs 
Cd by TT between A-f and A^. 

This paper is constructed as follows. First, in the next section we introduce the Hamiltonian and show the breaking 
of C/(l) symmetry by the interband scattering. In sec. |I| and |^ we investigate the effect of impurities on the 
bound and continuum states. It is shown that the bound states emerge only when ReA is positive and that the 
binding energy depends on the phase of the order parameter. We also show that the in the continuum states the 
energy gap is strongly suppressed when ReA = 0. In sec. M, we solve Bogoliubov-de Gennes equation for excitonic 
insulator numerically. The results support the analytic discussion above. We also find that, when is close to tt/2, 
magnetic moment of doped electrons is reduced. In the last section for summary and discussion, we discuss about the 
application of our results to the ferromagnetism in hexaborides. 
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II. BREAKING OF t/(l) SYMMETRY BY INTERBAND SCATTERING 

First we study the following hamiltonian 

k,(T k,k' .q.a,{T' k,k' ,q.a,{T' 



Here ak^a and bk^a are annihilation operators of electrons in the lower and upper band, e^ and e| are kinetic energy 
of electrons, and V and U are exchange and direct interaction between two bands. The nesting of Fermi surfaces 
between two the electron bands brings the formation of excitonic order parameteiQ. 
We introduce the interband scattering by impurity ^jj^irj as 

^imp = Y. ^imp 4a ha + c.c. (2) 

This interband scattering is important because this term breaks the f/(l) symmetry of the system. The hamilto- 
nian(P for the clean system is invariant under the U{\) x C/(l) transformation, a^ — > UiC^'^^ and bi — > 6^6*"^^. However, 
the interband scattering term reduces this symmetry to the diagonal U{1) symmetry because the U{1) x U{1) transfor- 
mation change the impurity term as V^jj^^aj^bicr -^ V[j^„al^bice^^'^^~'^'-\ unless (j)i ^ (j)2- This means that to discuss 

the impurity effect we must take the phase of excitonic order parameter, which corresponds to {a^b), into account. 

Here we note that in real materials this t/(l) symmetry is generally broken even without impurities due to other 
interactions which are not included in the Hamiltonianffl). The interband pair scattering term V a\^a\, _j3k' -abka 
is one of such interactions. Therefore we suppose the phase (/) is determined by such term. It is interesting to study 
the stable phase 4> and discuss the phase dynamics, however, this is beyond the scope of this paper. j.. 

Volkov et al. noticed that localized states with uncompensated spin appears in the doped triplet excitonic insulatom. 
They concluded that the bound states appears only when A is positive. Their result means that the phase of the 
order parameter determines the existence of bound state. However, they assume that A to be real, which is needed to 
be justified by further study. In this paper, we treat the order parameter as the complex variable. Here we introduce 
the excitonic mean field order parameter Ao-o-' and write down the mean-field hamiltonian in pure excitonic system 

HmY ^Y^>la(^ka + elbl^bka - Y^ {^aa'bl^,aka + ^l„,al^,bka) - tiY^al^'^ka +bl^bka)- (3) 

When we use the hamiltonian ([l|), order parameter Aa-a-/ is defined as 

^-' = ^E(4>fe^') -UY S.a'iaL^bka"). (4) 

k k,all 

This Hamiltonian is invariant under the transformation Oj — > aiC^'^^ , bi — > biC^'^'^ and Acra' ~^ e^^'^^~'^^\ It is convenient 
to assume the order parameter Aq-q-/ to be real by choosing the appropriate gauge. In this case, the impurity scattering 
term (||) shows phase dependence, 

^imp = E^impe^V>- + c.c. (5) 

In the following two sections, we use Hamiltonian (0) and (||). 

III. BOUND STATE FORMATION AT AN IMPURITY 

First, we study the bound state around the impurity in excitonic insulator. We here consider the one-impurity 
problem. For simplicity, we assume the order parameter Aa-a-' to be diagonal, Aa-a-' — AaS^a'- In the singlet excitonic 
insulator this assumption is correct. We note that in the triplet excitonic insulator we can justify this assumption if 
the three component of triplet order parameter A^ = X]q /3 ('''«) a/3 '^a/^ (* = 2;,2/,z) are degenerate. In this case, we 
can treat up-spin and down-spin separately. Then from (ph the Green's function in the mean-field theory is given by 

Assuming the impurity potential as point-like, we can calculate the Green's function Gaa'{k,u},k' .lo) using the 
T-matrix approximation, 

a" 



and 



After some calculation with the assumption e"(fc) = — e^(fc) and constant density of states, we get 

r.,, = F™P(i - J2 G"(fc)F"^^P)-\ (10) 



and 



>imp / A^e-"*' {-Lu - ii)e">' 



Y ^A2 - (cj + /i)2 V (-^ - ^^y ^ ^-e"^ ; ' ^ ^ 

here A^j? is the density of states at Fermi energy. 

The determinant of matrix D = \ — ^^, G"g.(fc)y^™P is given as 



^Al-{LO + y)^" ^Al-[u + llY' A2-(a> + /i)2 



= l + f^-2 /^-^"^^ , (12) 

where y = '"'^■F^imiD- From this expression, we find that when Re(V"Acos(/)) is positive D^^ has poles at w = 

— ^ ± W (1 i>2-|2 )^o- ^^'^ t'^o bound states appear. On the other hand, when Re(yAcos0) is negative, there 

exists no bound state. We conclude that the phase of the order parameter controls the energy and existence of bound 
states. If wc set A real and positive, there exist bound states at cj = —// ± j- — ^r^ A^. when = 0. Here we note that 
this bound state exists even if V-y^^„ is very weak. As (p changes from to 7r/2, the energies of bound states approach 
to oj = — /x ± A and at = 7r/2 the bound states are touch the edges of the continuum. 

It should be noticed there exists a great difference in the effect of impurities between the singlet excitonic insulator 
and the triplet one. In the singlet excitonic insulator, A| — Aj^ and both up-spin and down-spin bound states have 
the same energy. On the other hand, in the triplet excitonic insulator, A^ = — A| and the bound states of up-spin and 
down-spin electrons cannot exist at the same time. Therefore spins of the electrons doped into the triplet excitonic 
insulator align ferromagnetically. For (p = 7r/2, marginal bound states exist for both spin states and excitonic insulator 
remains paramagnetic. 

IV. REDUCTION OF ENERGY GAP BY IMPURITIES 

In the previous section, we studied the formation of a the bound state around an impurity. We found that a bound 
state appears only when Rc(V^Acos(/)) is positive. However, the impurity will also influence the continuum state. For 
example, it is well known that magnetic impurities in a superconductor reduce the energy gap. Zittartz discussed the 
effect of intraband scattering by impurities in excitonic insulator and found that it reduces the energy gap in exactly 
the same way as magnetic impurities do in a superconductorQ. In this section, we consider the reduction of energy 
gap by interband scattering term. I-. 

To study the suppression of the energy gap, we use the same technique as Abrikosov and Gor'kovB used for a 
superconductor. For simplicity we again assume that the order parameter is diagonal in spin space. We introduce the 
renormalized Green's function G 

G,,,{k,u;)^(^-f + f' , -^^ ) \„„,. (13) 

Here uj and A are the renormalized frequency and order parameter respectively, and we recovered the translational 
invariance of Green's function of continuum states by averaging over the impurity positions. 
Using the Born approximation, we get 



G,,,{k,iUy^ = G°,,{k,Uj)-^ - ^aa'{k,Uj) 



(14) 



and 



Eo-o-'(fc,Lj) = n 



imp 



d^k' 

(2^ 



l^aa" Lf cr"<T"' (.K , LJJ K^,,,^, 
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Integrating over fc' we obtain 



T,^^'{k,u) 



^imp^imp^^i- 



(i + A* 



-a) + /i 



-^aa' ; 



(15) 



(16) 



where we assumed e^ = — e^ and constant density of states as before. After some calculations we get the following 
equations for uj and A^- 



OJ = UJ + ■ 



w + M 



i^ + fiy 



(17) 



A, = A, - - 

T 



\A,\^-iu + fiy 



(18) 



Here r = l/rijj^^NpV-'^ . If = 0, these equations are the same as those for magnetic impurities in a superconductor. 

Therefore if (/> = 0, the energy gap is strongly suppressed by impurities. As (j) develops from 0, the energy gap increases 
and at = 7r/2, the energy gap becomes almost as same as that of pure excitonic insulator, because at = tt/2 the 
equations (17) and ( [l8|) have the same form as those for non-magnetic impurities in a superconductor. 
The density of states of electrons N{uj) is given by 



N{uj) 



i-Im|TrG(fc,c.)A^ 



_fk_ 
2tt""J "^"^'"'(2713) 
(i + /i 



= -iV^Im^(^ 



A,|2-(c:, + ^)2 



(19) 



For some typical values of (j), we plot the density of states in Fig. nl From this figure, wc find that the reduction of 
the energy gap is maximum at ^ = 0. This is consistent with the above discussion. 
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FIG. 1. phase dependence of DOS when r=100 



V. NUMERICAL SOLUTION OF BDG EQUATION 

From the analysis in the previous sections, we found that the phase of the order parameter determines the existence 
of the bound state as well as the reduction of the energy gap. These two effects are not independent. For example, 
for the finite concentration of impurities, the bound states will form an impurity band. On the other hand, if the 
energy gap of the continuum is strongly reduced, the energy of the bound state present for a single impurity model 
may be absorbed into the continuum for a finite impurity concentration. To study such possibilities, in this section 
we solve the Bogoliubov-de Gennes equations for excitonic insulator numerically. We confine ourselves to the triplet 
case, because in the triplet excitonic insulators impurity may invoke magnetic moment as discussed in section [11 . We 
also show that incomplete ferromagnetism can appear. 

We start from the following Hamiltonian 

%dG = J2 (-*"«l«J- - t'blb,,) - J2 {^^,<ra4a'a^a + C.C.) 
(i^j) ,(T i,a.a' 

+nmp E(^l°- + ^•^■) - ^- E «l«- ~^bY. ^l^- (20) 

i,(T i,a 

The first term is the kinetic term of electrons. Here we assume two-dimensional square lattice and the summation is 
carried out with nearest neighbor {i^j). In the theory of the excitonic insulator, it is often assumed that ef, = —t\. , 
so we take i° = —t^ = 1. In this case, e^ = — 2(cosfc^ + cosfcj^) — Ea and e\ = 2(cosfca; + cosfcj.) — E'f,. The second 
term is the interaction term. The order parameter Ai^g-cr' is determined from the Hamiltonian (nl) as 

Aw = V{alh„.) - U6„„. Y.'^al,Aa")- (21) 

a" 

While for [/ = triplet state and singlet state are degenerate, positive U suppresses singlet one. In the following 
numerical calculation we take U large enough to suppress the singlet state completely. The third term describes the 
impurity scattering and the sum is carried out on the impurity site. The last two terms are the energy gap term, 
which determines the overlap of two bands. In the following, we set Ea = —Ei,. Under this assumption, the Fermi 
surfaces of two electron bands coincide when we dope no electron. The size of the Fermi surfaces depends on Ea- For 
Ea < —4:t, e^ < e^ for any k, and the system is a band insulator. On increasing £'", the Fermi surface emerges near 
k = (7r,7r) and the system becomes a semi-metal. When Ea > ^t it becomes a band insulator again, because e\. < e^. 



We solve equation (l2d) and (21) numerically keeping the number of electrons constant. First, we diagonalize the 
equation(|2y) and calculate the eigenenergy and eigenstates. Next we take the eigenstates with small eigenvalues, up to 
the number of electrons. The order parameter is calculated by equation (Ell). We repeat these procedures iteratively 
until the solution converges. We change the phase (j) by choosing the initial values of A, and check the phase after 
convergence. 

We set parameters U = 1.0, V = 2.0, V-y^^ = 1.0 and Ea — —Ef, — 0. The calculation is carried out on 17 x 17 
lattices. 

First we consider the case when the impurities are introduced into the non-doped excitonic insulator. In Fig. we 
present the triplet order parameter A2 axis when one impurity is introduced, setting the initial condition = 0, 7r/4 
and 7r/2. The order parameter is reduced at the impurity site. We note that in Fig. 0(b) differs from 7r/4. This is 
due to the fact that after iteration deviates from the initial values. However, from Fig. g we find that the difference 
is small. 
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FIG. 2. Real and imaginary part of triplet order parameter when 



= 0,0 = 7r/4 and = tt/2. 



To see the reduction of energy gap, we need finite concentration of impurities. In Fig. ^ we plot the single particle 
spectrum near e = when 10 impurities are introduced periodically. We find that the bound states appear around 
E = 0.2 when (f> = 0. It should be noted that the all electrons in the bound state has up-spin, as shown in the right 
of the figure. The number of localized eigenstates with negative energy is the same as that of impurities. The energy 
of the bound states become large as (j) changes from to tt/2, and vanish at 7r/2, in consistent with our discussion. 
We also note that at (j) = 7r/4 the energy gap is larger than the case when (/> = 0. This is consistent with the result 
of section IV. We note that at = 7r/2 the gap seems to be smaller than the case when (j) = 7r/4. At i/i = 7r/2, the 



bound states touch the edge of the continuum and hybridize with the continuum. The energy of eigenstates is lowered 
by this hybridization. 
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FIG. 3. phase dependence of energy spectra near e = 0. phase (j) is set to 0, n/4 and n/2. 

spectrum of up-spin and down-spin electrons separately, when = 0. 



The right figure shows the 



Now wc turn to the case when electrons are doped. When doped into a triplet excitonic insulator, electrons will align 
ferromagnetically if ^ 7r/2. In the case of La2:Cai_2:B6, La acts as both the impurity and electron donor. Therefore 
we consider the case when the number of impurities and doped electrons are the same. In Fig.H we plot rii^i —rii^i when 



10 electrons and impurities are doped for (j> = 0, 77r/16, 7r/2. When (j) — we can see that locaUzed magnetic moment 
appears around the impurity. Total magnetic moment amounts to 10, which means that all doped electrons have spin 
up. On the other hand, when cj) — 7r/2 no magnetic moment appears. Therefore we can conclude when (f> = tt/2 the 
doped triplet excitonic insulator is paramagnetic. These results are also consistent with our discussion. In the middle 
of these two limits, there may exists incomplete ferromagnetism. In Fig.H(b), we show the same plot for (f> = Tvr/IG. 
We can see that magnetic moment become small, being 6 in total. This incomplete ferromagnetism is considered as 
follows. If the concentration of impurities is finite, associated impurity bound states form an impurity band. Since 
the bound state energy approaches A as (/> — > 7r/2, this impurity band has a finite overlap with the continuum state 
for (j) sufficient close to 7r/2. In this case ferromagnetic moment due to doped electrons will be reduced. 
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FIG. 4. Spatial pattern of magnetization when (a) 
the number of impurities are both 10. 



0, (b) — 77r/16 and (c) — n/2, when the number of electrons and 



VI. CONCLUSION AND DISCUSSION 



On the basis of the results of the preceding sections, we first discuss the weak ferromagnetism in doped hexaborides. 
As discussed in this paper, doped triplet excitonic insulator with impurities shows ferromagnetism if phase </> is between 
— 7r/2 and 7r/2. The incomplete ferromagnetism can also occur when the energy gap becomes smaller than the bound 
state energy. Here we discuss the possibility that the ferromagnetism of La-doped CaBg is caused by impurities. We 
assume that CaBg is the triplet excitonic insulator. When we dope La, La acts as impurity and cause interband 
scattering. As discussed in this paper, in the presence of interband scattering up-spin and down-spin does not 
compensate because of the emergence of localized moment if ^ 7r/2. Therefore ferromagnetism emerges, which 
become incomplete when the reduced energy gap is smaller than the highest energy of impurity band. Here we discuss 
the possibility of the impurity-induced ferromagnetism by this mechanism. 

First, we notice the ferromagnetism emerges only when excitonic order is spin-triplet. It is well known that direct 
interaction term suppresses the singlet order, while exchange interaction stabilize both singlet and triplet order by 
an equal amount. Therefore only triplet order parameter will be realized when direct interaction is large enough. It 
seems natural to assume that the exchange interaction is the same order as the direct interaction. As shown in sec.^, 
triplet excitonic insulator is realized when U and V are the same order. 

Second, the phase of order parameter must not be 7r/2, in other words, the minimum of the total energy is achieved 
when (j) ^ tt/2. In our discussion, we do not discuss the value of </> which minimize the total energy. In our numerical 
calculation the change of energy by phase is very small. As noted in sec.|| in real materials there exist some other 



interactions which is not discussed here, such as interband pair scattering. These interactions also depend on the 
phase. Because these interactions exist even without impurities, we think that the stable phase cf) will be determined 
by such term. Further study will be needed to determine the stable (j), however, there exists no reason to restrict 
(/) = 7r/2. 

Third, the minimum of the energy gap must be smaller than the energy of the impurity band. In our calculation 
it means that (p is close to tt/2. However, it should be noticed that there exists other origins of the reduction of the 
energy gap. For example, intraband scattering by impurities also decreases the energy gap, as discussed by ZittartzEl. 
In this paper we assume La acts as the interband scattering potential. However, there exists other impurities or 
disorder in La3;Cai_2;B6. These impurities and disorder will also acts as impurity potential. If these scattering is 
large enough, the energy gap will be reduced such that the energy of the bound states is larger than the gap. We 
conclude that the minimum of the gap can be smaller than the energy of the bound states. 

One may think that our model is too simplified to describe the La-doped hexaborides. In the following we briefly 
comment on the physics neglected in our model. 

First, we only consider one pair of Fermi surfaces, which shows completely nesting without doping. However, t 
band calculation of CaBg shows that there exist three pairs of imcompletely nested Fermi surfaces. Balents and Varm, 
discussed that in CaBg we must consider the possibility of the co-existence of gapped and ungapped Fermi surfaces. 
In this case, though our theory should be modified, ferromagnetism will occure. Our theory is based on the fact that 
Green's function depends on the phase. Even if gapless Fermi surfaces exist. Green's function Ga,(T'{k,uj,k' ,uj) will 
show the phase dependence. Therefore the occurence of ferromagnetism will be suggested in this case. On the other 
hand, the formation of LOFF state discussed by Barzykin and GorkovB gives severe damage to our theory. If LOFF 
state appears, order parameter A shows spatial modulation, such as A(r) = Aexp(zgr). In this case, total magnetic 
moment will diappear, because the direction of local magnetic moment will be random. 

Another important point we have neglected is the symmetry of the energy band. Band calculation shows that the 
conduction band and the valence band have different symmetries, Xg and X^. Though a simple substitution of Ga by 
La at a single site will reduce the symmetry, the rotational symmetry and parity around the impurity survive. This 
leads to the conclusion that interband scattering by La is forbidden by the symmetry. Therefore for our theory we 
must assume the symmetry is lowerd at La site. For example, if the impurity locates at less symmetric position, these 
symmetries break down and intraband scattering is permitted. 

We also note that the disorder in the lattice will also cause the intraband scattering. Though in the above discussion 
we have assumed that only La acts as interband scattering potential, interband scattering by disorder will also happen. 
In this case our theory does not modified, except that the number of impurities is different from that of doped electrons. 

The readers will also notice that intraband scattering is neglected in our model. However we can show that the 
bound state can survive even if intraband scattering exists. When one impurity causes both interband and intraband 
scattering, the energy of localized state is determined by 

(1 + ^i'mp - nmp)(^' - (^ + ^'f) + 2^imp(^ + /^) aV-^pAcos </. = (22) 

where V-' = TrNpV-' and V-' is the strength of intraband scattering. If \V-' \ < V^^-^p cos cj), equation ( p2[ ) has 
two bound states when Vjj^^p A cos > . In this case the discussion above is not qualitatively modified. On the other 



hand, \V-' | > Vjj^^ cos cj), equation (E2[) has one bound state, independent of the sign of VA cos (j). However here we 

note the energy of bound state differs between positive and negative A. Therefore we conclude the ferromagnetism 
will also survive, though some modification will be needed. 

So far there are no direct evidence of this impurity-induced incomplete ferromagnetism in the hexaborides. However 
we note that we can account for the experimental inconsistencies on the value of the saturation moment of La^jGai-ajBg 
by this scenario. For example, in the first paper of Fisk et al. the saturation moment is about 0.07^b per electron 
when X ~ 0.005, while Terashima et al. repcpts that the saturation moment is about 0.25 fj,B with same doping, 
more than three times larger than that of Fiskliy. We can explain this discrepancy by the difference of the density of 
impurities. More careful experiments on the magnetism will be needed to test the present theory. The measurement 
of the local magnetic moment will give much information about this. 

In summary, we have investigated the effect of the impurity in the excitonic insulator. The effect depends on the 
phase of the order parameter. The bound state appears when ReA is positive, and the energy of the bound state 
depends on the phase of the order parameter. On the other hand, the reduction of the energy gap of continuum state 
becomes maximum when ImA = 0. Numerical solution of BdG equation supports these discussion. Based on these 
results, we proposed a new mechanism of incomplete ferromagnetism and applied it to the hexaborides. Our theory 
qualitatively explains experiments, though further study will be needed. 
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